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Abstract 

Up to now, critical loads calculations for the Flemish Region were based on European background data of 
surrounding countries. A first attempt has been made to calculate and map critical loads for forest ecosystems in 
Flanders using available site-specific information. Values of current deposition were used to calculate and map 
exceedances. The lowest critical loads for acidification (697 eq ha -1 year -1 ) occur in the Campine and the north of 
Limburg where ecosystems largely consist of coniferous forests on poor sandy soils. The dominance of coniferous 
forest types in the Campine is also responsible for low critical load values for eutrophication (between 536 and 971 
eq ha -1 year -1 ). In 75% of the receptor points that have been considered an exceedance of the critical load for 
acidification is noted, primarily in areas with high S0 2 and NO x depositions, such as the north of the provinces East 
and West Flanders and Antwerp. The critical load for eutrophication is exceeded in all points considered. 
Exceedances are particularly high in coniferous forests in West Flanders, and in the north of the provinces of 
Antwerp and Limburg, where especially NH r depositions amount to high values. Data needed for the calculation of 
critical loads are still sparse in Flanders, e.g. for: (1) weathering rates of soil minerals; (2) interception and 
evaporation of forest ecosystems; and (3) uptake of N and basic cations by vegetation. This supplementary 
information will contribute to a further refining of the calculated critical loads, which constitute indispensable 
information in developing an emission abatement policy. © 2000 Elsevier Science B.V. All rights reserved. 
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1. Introduction 

Input of air pollutants such as S0 2 , NH X and 
NO x into forest ecosystems can result in damage 
to the vegetation both by direct effects on the 
forest canopy and by indirect soil-mediated 
processes of acidification and eutrophication in 
the soil. Direct harm to foliage by pollutants may 
contribute to the decline of forest vitality in Eu¬ 
rope. Besides these direct effects, the disturbance 
of the nutrient status of forest soils caused by soil 
acidification and nitrogen accumulation is of ma¬ 
jor importance (Roberts et al., 1989). 

Deposition of S0 2 , NH X and NO x can cause 
low pH values and high aluminium concentra¬ 
tions in the soil solution. Exchange of basic cations 
such as Ca 2+ , Mg 2+ , K + , and Na + by protons, 
adsorption of SO 2- on aluminium and iron hy¬ 
droxides with the release of OH~ into the solu¬ 
tion (Johnson, 1984) and acid-base reactions in 
the soil solution are fast processes which can 
temporarily counteract acidification. The release 
of basic cations through mineral weathering, on 
the other hand, is slow but provides a virtually 
inexhaustible source of buffering capacity (De 
Vries, 1994). Atmospheric deposition of basic 
cations is also a long-term source of buffering 
capacity. Besides their capacity to replace acidity 
from the sorption complex, basic cations may also 
improve the nutrient status of the ecosystem by 
balancing a high nitrogen input. 

If rates of acid deposition input exceed the 
rates of acid neutralising processes, the soil will 
acidify. Excessive concentrations of Al 3+ and H + 
in the rooting zone can then exert a toxic effect 
on the roots of trees and may inhibit basic cation 
uptake, causing unfavourable ratios of minerals in 
the plant tissue. 

Nutrient input into any ecosystem may result in 
eutrophication. Nitrogen is a growth-limiting fac¬ 
tor. Initially, an extra supply of N will induce an 
increase in ecosystem biomass production. Upon 
further accumulation of inorganic nitrogen, nitro¬ 
gen-tolerant species will gain competitive advan¬ 
tage over less tolerant species and the possibility 
of nitrogen leaching will increase. An additional 
effect is the increased susceptibility to frost and 


fungal diseases related to high leaf nitrogen con¬ 
tents (De Vries, 1988, 1990a,b). 

In the frame of the Geneva Convention on 
Long Range Transboundary Pollution, ecosys¬ 
tems in Europe are currently being mapped re¬ 
garding their sensitivity to acidification and eu¬ 
trophication in co-operation between members of 
the UN/ECE Task Force on Mapping (Downing 
et al., 1993). From physical, chemical and 
geochemical properties of the ecosytem, critical 
loads are calculated. The critical load is defined 
as The highest deposition of a compound that will 
not cause chemical changes leading to long-term 
harmful effects on ecosystem structure and func¬ 
tion’ (Downing et al., 1993). By comparison of the 
calculated critical loads with actual deposition 
rates, so-called areas of exceedance can be lo¬ 
cated. 

The critical loads concept is an important in¬ 
strument in assessing the potential for acidifica¬ 
tion stress imposed on forest soils by anthro¬ 
pogenic depositions. It is also an indispensable 
tool to evaluate the efficiency of proposed emis¬ 
sion-reduction strategies (Posch et al., 1993) The 
purpose of determining critical loads is to set the 
goal for future deposition rates in order to pro¬ 
tect the environment. The objective of this study 
was to perform an initial trial to derive critical 
loads, using specific information for Flanders 
(Belgium). 

2. Materials and methods 

2.1. Computation 

To ensure compatibility between calculated val¬ 
ues from different National Focal Centres in Eu¬ 
rope, the ‘Co-ordination Centre for Effects (CCE)’ 
hosted by the RIVM (National Institute of Public 
Health and Environmental Protection, The 
Netherlands) proposed a standard protocol for 
the calculation and mapping of critical loads. The 
detailed procedure is described in the CCE-Status 
Report 1993 (Posch et al., 1993). The choice of 
the computational approach depends on the input 
data requirements, the complexity of use and the 
uncertainty tolerated on the calculated value. In 
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this study, the steady state mass balance model 
was used to derive the critical loads for forest 
soils only, as the procedure assumes that soils are 
subjected solely to natural acidification processes. 
The critical load for acidity, CL(A), considers 
acidifying effects of sulfur and nitrogen deposi¬ 
tion. Critical load for nutrient nitrogen, Cl nut (N), 
accounts for both the acidifying and the eutrophy- 
ing effects of atmospheric deposition of nitrogen. 
Exceedances of these critical loads are calculated 
as the difference between the actual loads and 
the critical loads. 

2.2. Data acquisition 

2.2.1. Receptor points 

The critical load is related to a receptor, i.e. the 
ecosystem that is subjected to the harmful effects 
of atmospheric deposition. Selection of these re¬ 
ceptors is the basis of the critical load calcula¬ 
tions (Hootsmans and Van Uffelen, 1991). Recep¬ 
tor points were located on a grid of 5 km X 5 km 
after an overlay of The Flemish Forest Inventory 
Map (FFIM) on the Soil Association Map of 
Flanders (Table 1). 

The lowest critical load and the highest ex¬ 
ceedance within every grid cell was then at¬ 
tributed to the whole cell. By this approach, every 
ecosystem within a grid cell is protected by this 
critical load and the worst exceedance situation 
within every cell is displayed. 

2.2.2. Weathering rates 

Because more specific data on the release of 
basic cations through mineral weathering are not 
yet available in Flanders, weathering rates were 
assessed according to the method proposed by 
Hettelingh and De Vries (1992). Soil types are 
classified into four parent material classes and six 
texture classes. A combination of both classes 
indicates the weathering class to which the recep¬ 
tor point belongs. To each weathering class, base 
cation release rates were assigned on the basis of 
mineralogy, according to Sverdrup and Warfvinge 
(1988). Weathering rates were calculated for a 
root zone of 50 cm. This is the rooting depth for 
trees, commonly assumed in critical load compu¬ 
tation. 


Table 1 

Soil and forest types considered and the proportion of the 
resulting forest-soil combinations (receptor types) in the total 
number of receptor points 


Soil type a 

Lithology 

Re Calcaric Regosol 

Recent dune sand 

Jc Calcaric Fluvisol 

Recent alluvium 

Lg Gleyic Luvisol 

Loess 

Lo Orthic Luvisol 

Loess, tertiary sand and clay 

Dd Dystric 

Podzoluvisol 

Sandy loess, tertiary green sand 

Ph Humic Podzol 

Pleistocene aeolian sand 

Po Orthic Podzol 

Pleistocene aeolian sand, 
tertiary sand 

Forest type 

Description 

1 Deciduous forest 

< 20% coniferous trees 

2 Mixed deciduous 

20-50% coniferous trees 

3 Mixed coniferous 

20-50% deciduous trees 

4 Coniferous 

< 20% deciduous trees 


Receptor type 

(%) 

Receptor type 

(%) 

Lo-1 

29.1 

Po-1 

0.8 

Ph-4 

24.2 

Lg-1 

0.6 

Ph-1 

14.9 

Po-3 

0.5 

Dd-1 

10.1 

Dd-2 

0.3 

Ph-3 

6.1 

Lo-3 

0.3 

Po-4 

4.1 

Rc-1 

0.3 

Ph-2 

2.5 

Jc-1 

0.2 

Dd-4 

2.0 

Lg-3 

0.2 

Lo-4 

2.0 

Lg-4 

0.2 

Dd-3 

0.8 

Po-2 

0.2 

Lo-2 

0.8 




a FAO-UNESCO (1981). 


2.2.3. Precipitation 

Mean annual precipitation over 10 years 
(1986-1995) was derived from data of five weather 
observation stations in Flanders. Each receptor 
point was assigned the value corresponding to the 
nearest observation station. 

2.2.4. Interception 

The fraction of the precipitation which is inter¬ 
cepted by coniferous and deciduous trees differs 
due to the difference in leaf carrying period and 
leaf type. Values were derived from Hootsmans 
and Van Uffelen (1991). 
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Table 2 

Minimum and maximum total depositions (S0 2 + NO x + NH X ) 
and N depositions (NO x +NH X ) per grid cell 



Total dep. 

Total N d 


(eq ha -1 year -1 ) 

(eq ha year : ) 

Min 

2780 

1460 

Max 

4110 

2320 


2.2.5. Evapotranspiration 
Evapotransporation in forest stands is difficult 

to assess as it is influenced by many factors such 
as leaf area, leaf carrying period, soil moisture 
content, etc. In Flanders, annual evapotranspira¬ 
tion ranges between 270 and 370 mm year 1 (Van 
Miegroet and Dua, 1985). A value of 320 mm 
year 1 was therefore assumed. 

2.2.6. Deposition 

Total deposition per component was calculated 
for a 5 km X 5 km grid. Filtering factors, that 
account for the increased capture by forest filter¬ 
ing, were taken from Erisman (1990) for conifer¬ 
ous trees and from Ivens et al. (1988) for decidu¬ 
ous trees. 

Atmospheric deposition of basic cations may be 


of importance in soils where unweathered min¬ 
eral contents and hence release of basic cations is 
small. These data are indispensable for the calcu¬ 
lation of critical loads. Unfortunately, data on dry 
and wet basic cation deposition are as yet not 
available in Flanders. Total deposition of basic 
cations can be approximated by the total Ca 2+ 
deposition (De Vries, 1994). On the assumption 
that Flemish conditions are similar to those in 
The Netherlands, 1-year measurements of the 
mean total Ca 2+ deposition near 10 forest plots 
(Van Den Berghe et al., 1991) were used to 
estimate deposition of basic cations. 

2.2.7. Uptake by vegetation 

The mean annual uptake of nitrogen and basic 
cations were estimated from literature data (De 
Vries, 1990a,b). The uptake is strongly related to 
tree species, forest management and site quality. 
For coniferous and deciduous wood, mean nitro¬ 
gen uptakes were estimated to be 300 and 500 eq 
ha -1 year -1 , respectively. Hootsmans and Van 
Uffelen (1991) observed that the uptake of nitro¬ 
gen and basic cations differs more between forest 
types than between soil types. 
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Fig. 1. Deposition of S0 2 . 
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Fig. 2. Deposition of NO^ 


3. Results and discussion 

3.1. Deposition 

Deposition patterns are remarkably different 


between components (Figs. 1-3). The dominant 
source of S0 2 and NO x is the combustion of 
fossil fuels. Values are high in urbanised and 
industrialised areas such as the area between 
Ghent, Antwerp and Brussels. Intensive animal 
husbandry, on the other hand, is an important 
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Fig. 3. Deposition of NH X . 
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Table 3 

Minimum and maximum total depositions (S0 2 + NO x + NH X ) and N depositions (NO x + NH X ) per ecosystem 



Total dep. 

1 1 
(eq ha - year - ) 

Type 

1 1 

(eqha - year - ) 

Total N dep 

Type 

Min 

2613 

Rc-1 

1253 

Rc-1 

Max 

6228 

Po-3 

2568 

Lo-3 


Table 4 


Minimum critical load for acidification and nutrient nitrogen per receptor type (eq ha 1 year ] ) 

Soil 

Forest type 







1 


2 


3 


4 


CL(A) 

CL nut (N) 

CL(A) 

CL nut (N) 

CL(A) 

CL nut (N) 

CL(A) 

CL nut (N) 

Jc 

6053 

874 







Re 

6000 

812 







Dd 

1512 

804 

1573 

808 

1248 

538 

1248 

538 

Lg 

1976 

819 



1627 

573 

1627 

544 

Lo 

1512 

804 

1678 

808 

1377 

538 

1336 

536 

Ph 

1012 

804 

1012 

804 

697 

536 

697 

536 

Po 

1071 

816 

1226 

819 

836 

542 

697 

536 


source of NH X . Indeed, high deposition rates are 
encountered in the West of Flanders. Total depo¬ 
sition values on the grid cells ranged from 2780 to 
4110 eq ha -1 year -1 (Table 2). 


Mean values for total deposition per ecosystem 
include filter factors to account for the effects of 
the forest leaves in retaining deposition (Posch et 
al., 1993). Values were lowest for combinations of 



CL(A) (eq/ha.yr) 


0 25 50 Kilometer 





500-1000 

1000-1500 

1500-2000 

>2000 

No receptor points 


Fig. 4. Critical loads for acidification [CL(A)]. 
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N 



ES3 1000-1500 

M 1500-2000 

■i >2000 

I I No receptor points 

Fig. 5. Exceedances of the critical load for acidification [E(A)]. 


calcaric Regisols with deciduous forest. Mixed 
coniferous forests on orthic Podzols were subject 
to the highest mean total depositions (Table 3). 

Total N deposition was lowest for coniferous 
forests on calcaric Regisols. The highest mean 
total N depositions were calculated for orthic 
luvisols under mixed coniferous forest. 


3.2. Critical loads 

The spatial distributions of the critical loads for 
acidification and exceedances over the Flemish 
Region are shown in Figs. 4 and 5. Figs. 6 and 7 
show the spatial distribution of the critical loads 



50 Kilometer 


CL(N) (eq/ha.yr) 

0-600 
BB 600-900 
>900 

1 1 No receptor points 


Fig. 6. Critical loads for nutrient nitrogen [CL nut (N)]. 
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and the exceedances of the critical loads for 
nutrient N, respectively. The blank grid cells on 
the maps contain no receptor points. In Table 4 
the lowest values of the critical load for acidifica¬ 
tion [CL(A)] and for nitrogen per receptor type 
are given. 

The highest values of the critical loads for 
acidification (6053 eq ha -1 year -1 ) were found 
for calcareous soils under deciduous forest (soil 
types Jc and Re, forest type 1). The estimated 
rate of base cation release from mineral weather¬ 
ing is high in these soils and therefore provides a 
long-term buffering capacity against acidification. 
Most vulnerable are combinations of coniferous 
forests and poor sandy soils with a low long-term 
reserve of neutralising base cations (soil types Ph 
and Po, forest types 3 and 4). The lowest critical 
loads were determined for The Campine and the 
north of Limburg where the ecosystems largely 
consist of very sensitive coniferous forests on 
poor sandy soils. 

The critical loads for nitrification [Cl nut (N)] (Fig. 
6) were calculated according to the constant deni¬ 
trification method (Posch et al., 1993). The higher 
sensitivity of conifers to euthrophication is due to 
the nitrogen uptake, which is lower than for de¬ 
ciduous tree species. The dominance of conifer¬ 


ous forest types in The Campine is responsible 
for the low Cl nut (N) in this region. Values per 
receptor type ranged from 536 to 874 eq ha -1 
year -1 (Table 4). 

3.3. Exceedances 

Critical loads for acidification were exceeded 
for 75% of the receptor points. The highest ex¬ 
ceedances of CL(A), Ex(A) amounted to 6630 eq 
ha -1 year -1 in the most sensitive areas. These 
are mainly found in the north-west of the Flemish 
Region (Fig. 5). Highest exceedances were calcu¬ 
lated for coniferous forests on poor sandy soils 
situated in areas with high S0 2 and NO^ deposi¬ 
tion such as the north of the provinces East and 
West Flanders and Antwerp. One reason for their 
larger sensitivity is that forest filtering factors are 
significantly lower for conifers (1) than for decid¬ 
uous forest (1.8). NO x deposition contributes less 
to the exceedance as the filtering factor of both 
forest types for this compound is rather low. 

Values for exceedances of nitrification varied 
between 153 and 3501 eq ha -1 year -1 . The 
amount of nitrogen disposed and the forest type, 
in terms of filtering factors, is of great influence. 
‘Black points’, with exceedances above 2000 eq 



0 25 50 Kilometer 



Enut(N) (eq/ha.yr) 



0-500 

500-1000 

1000-1500 

1500-2000 

>2000 
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Fig. 7. Exceedances of the critical load for nutrient N[E nut (N)]. 
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ha -1 year -1 , are situated mainly in conifer stands 
in West Flanders, the north of the provinces 
Limburg and Antwerp, where high depositions of 
NH X were measured (Fig. 3). Again, the influence 
of NO a . deposition is rather limited due to the low 
forest filtering factors of trees for this compound. 

3.4. Uncertainty 

Uncertainties on critical loads and their ex¬ 
ceedances are caused by the application criteria, 
the available input data and the method applied. 
Due to the large-scale approach, the criteria ap¬ 
plied in the method are strongly generalised. To 
enable the comparison of critical loads on a Eu¬ 
ropean scale, it was agreed to use the same crite¬ 
ria in national critical load calculations of the 
national focal centers. The uncertainty on the 
input data is not negligible. All data originated 
from literature reviews and research performed in 
the past. If specific data were missing, values from 
neighbouring countries were used. In the model, 
several simplifications and assumptions were 
made. The accuracy of the model depends upon 
the extent to which these are justified. 

4. Conclusion 

The lowest critical loads for acidification (697 
eq ha -1 year -1 ) were calculated for the Campine 
and the north of Limburg where ecosystems 
largely consist of very sensitive coniferous forests 
on low-fertile sandy soils. The dominance of 
coniferous forest types in the Campine is also 
responsible for low critical load values for eu¬ 
trophication (between 536 and 971 eq ha -1 
year -1 ). In 75% of the receptor points considered 
an exceedance of the critical load for acidification 
is noted, primarily in areas with high S0 2 and 
NO^ depositions, which is the case in the north of 
the provinces, East and West Flanders and 
Antwerp. The critical load for eutrophication is 
exceeded in all points considered. ‘Black points’ 
are situated in coniferous forests in West Flan¬ 
ders, north Limburg and in the north of the 
province of Antwerp, where especially NH X depo¬ 
sitions amount to high values. 


The calculated values indicate trends and pro¬ 
vide a good initial estimation of the spatial vari¬ 
ability of ecosystem sensitivity to acidification and 
eutrophication in Flanders. Since political deci¬ 
sions concerning deposition abatement strategies, 
on a national as well as on a European scale, are 
based on the accurate regional assessment of 
critical loads, additional research will be needed 
to complete available databases and to update 
computations. In particular, data on (1) mineral 
weathering rates; (2) interception of rainfall; (3) 
evapotranspiration; and (4) nutrient uptake by 
forest vegetation in particular need to be refined 
and expanded in the future. 
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